Line Drawing with the 
NS32CG16;NS32CG16 
Graphics Note 5 

1.0 INTRODUCTION 

The Bresenham algorithm, as described in the "Series 
32000® Graphics Note 5" is a common integer algorithm 
used in many graphics systems for line drawing. However, 
special instructions of the NS32CG16 processor allow it to 
take advantage of another faster integer algorithm. This ap- 
plication note describes the algorithm and shows an imple- 
mentation on the NS32CG16 processor using the SBITS 
(Set BIT String) and SBITPS (Set BIT Perpendicular String) 

instructions. Timing for the DRAW LINE algorithm is given 

in Tables A, B and C of the Timing Appendix. The timing 
from the original Bresenham iterative method using the 
NS32CG16 is given in Table D. 

The bit map memory conventions followed in this note are 
the same as those given in the NS32CG16 Reference Man- 
ual and Datasheet, and all lines drawn are monochrome. 
Series 32000 Graphics Note 5, AN-524, is recommended 
reading. 

2.0 DESCRIPTION 

All rasterized lines are formed by sequences of line "slices" 
which are separated by a unit shift diagonal to the direction 
of these slices. For example, the line shown in Figure 1 is 
composed of 7 slices, each slice separated by a unit diago- 
nal shift in the positive direction. Notice that the slices of the 
line vary in length. The algorithm presented in this note de- 
termines the length of each slice, given the slope and the 
endpoints of the line. 

Depending on the slope of the line, these slices will extend 
along the horizontal axis, the vertical axis or the diagonal 
axis with respect to the image plane (i.e., a printed page or 
CRT screen). If the data memory is aligned with the image 
plane so that a positive one unit horizontal (x-axis) move in 
the image plane corresponds to a one bit move within a byte 
in the data memory, and so that a positive one unit vertical 
(y-axis) move in the image plane corresponds to a positive 
one "warp" (warp = the number pixels along the major axis 
of the bit map) move within the data memory, then the 
SBITS and SBITPS instructions can be used to quickly set 
bits within data memory to form the line slices on the image 
plane, as explained in section 3.1. For long horizontal lines, 
the MOVMP (MOVe Multiple Pattern) instruction is more ef- 
ficient than SBITS. This instruction is discussed in section 

3.1 and in the NS32CG16 Reference Manual. 

2.1 Derivation of the Bresenham SLICE Algorithm 

For the moment, consider only those lines in the X-Y coordi- 
nate system starting at the origin (0,0), finishing at an inte- 
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ger end point (x,y) and lying in the first partial octant, as in 
Figure 2. (The analysis will be extended for all lines in sec- 
tion 2.2.) The equation for one such line ending at (A,B) is: 

y = mx, 
where 

m = B/A 
is the slope of the line. Note that because the line lies in the 
first partial octant, A > 2B > 1 . 
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FIGURE 2 

Each pixel plotted can be thought of as a unit square area 
on a Real plane (Figure 3). Assume each pixel square is 
situated so that the center of the square is the integer ad- 
dress of the pixel, and each pixel address is one unit away 
from its neighbor. Then let A, represent the X-coordinate of 
the pixel, as shown in Figure 3. The value of Y at A, is: 

y = (B/A)Ai 
where y is Real. 

Since the address of each pixel plotted must have corre- 
sponding integer coordinates, the closest integer to y is ei- 
ther the upper bound of y or the lower bound. (Recall that 
upper and lower bounds refer to the smallest integer greater 
than or equal to y and the largest integer less than or equal 
to y respectively.) The original Bresenham algorithm was 
based on this concept, and had a decision variable within 
the main loop of the algorithm to decide whether the next 
y\ + -\ was the previous y, (lower bound) or y, -I- 1 (upper 
bound). For the SLICE algorithm, we are only concerned 
with when the value changes to yi + 1 , and the length of 
the previous slice up to that point. 
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The line from (0,0) to (45,6) is a first octant line with run lengths 3-7-6-7-6-7-3. Notice that a pixel is plotted before the run begins so that the actual number of 
pixels plotted is equivalent to the run length -HI. 

FIGURE 1 

Series 32000® is a registered trademarl< of National Semiconductor Corporation. 
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Y is incremented when the location of the half point Is beyond Aj, or when 
the true value of Y at Aj + i Is greater than Yj + Va. 

FIGURES 

In order for y, to be incremented along the Y-axis, the true 
value of real y at A, + 1 must be greater than or equal to the 
halfway point between y, and yi + 1 (Figure 3). If we let i 
increment along the Y-axis, then this half point occurs 
when: 

y= 1/2 + yi 
Or, because y, = i when incrementing along the Y-axis, 

y = (1 + 2i)/2. 
The real value of x at this point is: 

X = A(1 + 2i)/2B 
using X = (1 /m)y. The lower bound of this value of x repre- 
sents the x-coordinate of the pixel square containing the 
half point. 

Letting A, and Aj + i be two integer values of x where the 
real value of y is greater than or equal to the half point value 
yi + 1 /2 (Figure 4), then the run length extends from (A, + 
1,1 + 1) to (A| + i, I + 1). The run length can then be 
calculated as: 

Hi + 1 = Ai + 1 -Ai-1 

for i = 0,1 (B-2). Using the equation for x above, we 

can now better define A, as: 

Aj = (A/2B) + (lA/B). 

This equation has two real-valued divisions which are not 

suitable for an Integer algorithm. However, the equation can 

be broken down so that it only involves an integer-valued 

division and its integer remainder, which Is more efficient for 

processing. To do this we must define some intermediary 

integer values: 

Q = lower[A/B] I Lower bound of inverted slope! 

R = gl* (Integer residue of A modulo B I 

M = lower[A/2B] (Can also be defined as Q/2) 

N = 2bI* (Integer residue of a modulo 2B! 

T| = 2Bl(N + 2iR) ( Integer residue of (N + 2iR) 

modulo 2B! 
Note: aIs = B + A •lower[A/B]. 




Using the above values we can now define A, as, 

A; = (M + N/2B) + (IQ + IR/B) 

Ai = M + iQ + (N + 2iR)/2B 
Therefore, substituting A; and Aj + i into the equation for 
Hj + i, the intermediate horizontal lengths are, 

Hi + 1 = Ai + 1 -Ai-1 

Hi + 1 = (M + (l + 1)Q + lower[(N + 2(i + 1)R)/2B]| - 
(M + IQ + lower[(N + 2IR)/2B]) - 1 
Hi + 1 = Q + lower[(N + 2iR)/2B + 2R/2B] - 
lower[(N + 2iR)/2B] - 1 
Hi + 1 = Q - 1 + lower[(Ti + 2R)/2B] 
Analyzing the term lower[(Ti + 2R)/2B] It is shown that If T, 
+ 2R > 2B then the term becomes 1 , otherwise it becomes 
0. This is due to the definition of residue and modulo. The 
term T, is defined as: 

(N + 2iR) - 2B(lower[(N + 2IR)/2B]), 
which means that < T, < 2B. The same Is true for R: 

R = A - B(lower[A/B]), 
so that < 2R < 2B. Therefore, 

< Ti + 2R < 4B 
and, 

< (Ti + 2R)/2B < 2. 
The only possible integer values for this term are and 1. 
The term will equal if Ti + 2R < 2B, and It will equal 1 
when Ti + 2R > 2B, and H| + i will equal Q. The decision 
variable can now be defined as 

testvar = T, + 2R-2B. 
If testvar > then the horizontal run length is Q; If testvar < 
then the run length is Q-1. 

Looking again at the definition of Tj, a recursive relationship 
for the testvar can be formed. 
Tj + i = (N + 2R(I + 1)) - 2B(lower[(N + 2R(i + 1))/2B] 
Ti + 1 = (N + 2iR + 2R) - 2B(lower[(N + 2iR + 
2R)/2B] 
Since, as shown above, < (Ti + 2R)/2B < 2 then low- 
er[(Ti + 2R)/2B] < 1. In fact, if Ti + 2R < 2B then low- 
er[(Ti + 2R)/2B] = 0, and If T, + 2R > 2B then lower[(Ti 
+ 2R)/2B] = 1. Therefore, letting Tq = N, 
Ti + 1 = Ti + 2R if (Tj + 2R) < 2B 
Ti + 1 = Tj + 2R-2B if (Ti + 2R) > 2B. 
This gives the recursive relationship for testvar: 
testvar 1 + 1 = testvar, + 2R 
Hi = Q - 1 
if testvar , < 0. And, If testvar , > 0: 

testvar i + 1 = testvar i + 2R-2B 
Hi = Q. 
These recursive equations allow the intermediate run 
lengths to be easily calculated using only a few additions 
and compare-and-branches. 
The initial run length is calculated as follows: 

Ho = Ao = lower[A/2B] = M + lower[N/2B] = M. 
The final run length is similarly calculated as: 
Hf = M - 1 if N = else Hf = M. 



Run length Is calculated as An 
Is 1. 
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FIGURE 4 



Thus, the SLICE algorithm calculates the horizontal run 
lengths of a line using various parameters based on the first 
partial octant abscissa and ordinate of the line. The algo- 
rithm is efficient because it need only execute its main loop 
B times, which is a maximum of A/2, if A is normalized for 
the first partial octant. Compare this with the original Bre- 
senham algorithm which always executes its main loop A 
times. 

2.2 Extended Analysis for All Other Lines 

In section 2.1 the SLICE algorithm was derived for lines 
starting at the origin and contained within the first octant (B 
< 2A). The algorithm is easily extended to encompass lines 
in all octants starting and ending at any integer coordinates 
within the pre-defined bit map. The only modifications nec- 
essary for this extension are those relating to the direction 
of movement and in defining the coordinates A and B. 
In order to extend the algorithm to cover all classes of lines, 
the key parameters used by the algorithm must be normal- 
ized to the first partial octant. Those parameters are the 
abscissa and ordinate displacements and the movement of 
the bit pointer along the line. The abscissa and ordinate 
displacements of the line are normalized to the first octant 
by calculating: 

delta x = Xf - Xg and delta y = yf - ys 
which represent the abscissa (delta x) and ordinate (delta y) 
displacements of the original line. Then, the first octant 
equivalents of A and B will be: 

A = maximum (|deltax|,|deltay| I 

B' = minimum (|deltax|,|deltay| I 

B = minimum (B', A — B'l 
The next step In normalizing the line for the first octant Is to 
assign the correct value to the movement parameters. A 
line In the first octant and starting at the origin always has 
horizontal run lengths in the positive direction along the X 
(major) axis, and has diagonal movement one unit in the 
positive X direction and one unit In the positive Y (minor) 
direction. Since the SLICE algorithm calculates the run 
lengths independent of direction, variables can easily be de- 
fined which contain the direction of movement for each slice 
and each diagonal step within the different octants. 
Lines of different angles starting at the origin have slices of 
different angles. For example, a line of angle between 22.5 
degrees and 45 degrees has run lengths that are diagonal, 
not horizontal, and the direction of the diagonal step is hori- 
zontal, not diagonal. Because of this characteristic, it is con- 
venient to break the 8 octants of the X-Y coordinate system 
into 16 sections, representing all of the partial octants. 
Then, re-number these partial octants so that they form new 
octants as In Figure 5. These redefined octants represent 
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split. The origin is at the center of the drawing. Setting DELX positive on all 
lines makes opposite octants equivalent in the table below. 
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each of the eight angle classes of lines. For example, the 
lines In octants 3 and 7 are composed of diagonal (45 de- 
gree) slices in either the positive or negative direction, and 
have diagonal step in the vertical position. Lines In octants 4 
and 8 have run length slices in the vertical direction with 
diagonal steps in the horizontal direction with respect to the 
X-Y plane. 

In conclusion, the SLICE algorithm calculates successive 
run lengths in the same manner for lines in each octant. The 
only difference between the octants Is the direction of 
movement of the bit pointer after each successive run 
length is calculated. The run lengths and diagonal steps for 
each octant are given in Table I. Figure 5 shows the octants 
used by the SLICE algorithm. 

3.0 IMPLEMENTATION OF SLICE USING 
SBITS, SBITPS AND MOVMP 

The NS32CG16 features several powerful graphics instruc- 
tions. The SLICE algorithm described by this application 
note is Implemented with three of these instructions: SBITS, 
SBITPS and MOVMP. The SBITS Instruction allows a hori- 
zontal string of bits to be set, while the SBITPS Instruction 
can set vertical or diagonal strings of bits. The MOVMP in- 
struction, not detailed in this application note, can be used 
to set long strings of bits faster than SBITS when the length 
is more than 200 bits in the horizontal direction. The 
BIGSET.S routine given in the appendix uses this Instruction 
in conjunction with SBITS for long lines. These are very use- 
ful instructions for the SLICE run length algorithm, as will be 
shown in section 3.2. 



TABLE I 



OCTANT 


DELA 


DELB 


DIAGONAL MOVE 


RUN LENGTH 


1 &5 


DELX 


|DELY| 


1 + (±WARP) 


+ HORZ 


2&6 


DELX 


DELA-|DELY| 


+ 1 


±DIAG 


3&7 


|DELY| 


DELA-DELX 


±WARP 


±DIAG 


4&8 


|DELY| 


DELX 


+ 1 


+ WARP 



If DELX < then the starting and ending coordinates are swapped. This simplifies initialization. 



3.1 SBITS and SBITPS Tutorial 
SBITS: 

SBITS (Set BIT String) sets a string of bits along the hori- 
zontal axis of a pre-defined bit map. The instruction sets a 
string of up to 25 bits in a single execution using four argu- 
ments pre-stored in registers RO through R3. 

RO = (32 bits) Base address of bit-string destination. 

R1 = (32 bits, signed) Starting bit-offset from RO. 

R2 = (32 bits, unsigned) Run length of the line segment. 

R3 = (32 bits) Address of the string look-up table. 
The value of the bit offset is used to calculate the bit num- 
ber within the byte, assuming that the first bit is bit and the 
last bit Is bit 7. A maximum of 7 for the starting bit number 
added to a maximum of 25 for the run length requires a total 
of 32 bits. SBITS calculates the destination address of the 
first byte of the 32-blt double word to contain the string of 
set bits by the following: 

Destination Byte = Base Address -I- Offset DIV 8. 
Then, the starting bit number within the destination byte is: 

Starting Bit = Offset MOD 8. 
SBITS instruction then calculates the address for the 32-bit 
double word within the string look-up table (found in the 
NS32CG 1 6 manual) which will be OR'ed with the 32-bit dou- 
ble word whose starting byte address Is Destination Byte, as 
calculated above. The table is stored as eight contiguous 
sections, each containing 32 32-bit double words. Each of 
the eight sections corresponds to a different value of Start- 
ing Bit (Offset MOD 8), which has a possible range of 
through 7. The 32 double words in each section correspond 
to each value of the run length (up to 25) added to the 
starting bit offset. 
example: 

Register Contents 
before after 



RO = 1000 
R1 = 235 
R2 = 16 
R3 = $stab 



RO = 1000 
R1 = 235 
R2 = 16 
R3 = $stab 

1029 



Destination Address = 1000 + (235 DIV 8) ^ 

Starting Bit = 235 MOD 8 = 3 

Table Address = Sstab + 4*(16 + (32*3)) = $stab + 448 

bytes 
32-bit Mask = 0x0007FFF8 

This mask value Is OR'ed with the 32-bit double word start- 
ing at byte address 1029 decimal. Notice that the mask 
0x0007FFF8 leaves the first 3 bits and the last 13 bits 
alone. Thus, a string of 16 bits Is set starting at bit number 3 
at address 1029 decimal. The contents of the registers are 
unaffected by the execution of the SBITS Instruction. 
Since the SBITS instruction can set up to 25 bits in one 
execution, the run length In R2 can be compared to 25, and 
a special subroutine executed If it exceeds 25 bits. The sub- 
routine will set the first 25 bits, then subtract 25 from the run 
length, and compare this to 25 again. This process Is re- 
peated until the run length is less than 25, in which case 



the remaining bits are set and the subroutine returns. The 

DRAW LINE algorithm implemented In this application 

note uses this method for strings of bits to be set less than 
200. For horizontal lines greater than 200 pixels in length, 
the BIGSET routine is more efficient, as described below. 
BIGSET: 

The utility program BIGSET.S is used to draw longer lines, 
more than 200 pixels in length, more efficiently than SBITS. 
BIGSET.S, which is given in the appendix, uses the MOVMP 
instruction (MOVe Multiple Pattern) to set long strings of 
bits. Since MOVMP operates on double-word aligned ad- 
dresses most efficiently, the string is broken up into a start- 
ing string within the first byte, a series of bytes to be set, and 
an ending string which is the leftover bits to be set within the 
final byte. The starting and ending strings of bits, if any, are 
set using the SBITS table with an OR Instruction. 
SBITPS: 

SBITPS (Set BIT Perpendicular String) handles both vertical 
lines and diagonal lines. This instruction also requires four 
arguments pre-stored in RO through R3. RO, R1 and R2 are 
the Base Address, Starting Bit Offset and Run Length re- 
spectively, as for SBITS. R3, however, contains the destina- 
tion warp. 

Note: The Destination warp is the number of bits aiong the horizontai iength 
of the bit map, or the number of bits between scan iines. it is aiso 
referred to as the "pitch" of the bit map. Thus, a verticai one-unit 
move in the positive direction wouid require adding the vaiue of the 
warp to the bit pointer. A diagonai or 45 degree iine is drawn when the 
warp is incremented or decremented by one. 

The run length is a 32 bit unsigned magnitude. 

example: 

(Assume that the bit map is a 904 x 904 pixel grid.) 



before 



Register Contents 
after 



RO = 1000 RO = 1000 

R1 = 235 R1 = 235 + (150*904) = 135,835 

R2 = 150 R2 = 

R3 = +904 R3 = +904 

Destination Address = 1029 
Starting Bit Number = 3 
Run Length = 150 
Warp = + 904 

As In the example for SBITS, the Destination Address is 
1029, with Starting Bit Number = 3. Since the warp In this 
example is + 904 and the bit map is 904 x 904 bits, the line 
is vertical, has a length of 150 pixels and starts at bit num- 
ber 3 within the byte whose address is 1029 decimal. Unlike 
the SBITS instruction, the SBITPS alters registers R1 and 
R2 during execution. R1 is set to the position of the last bit 
set plus the warp. However, this is convenient for drawing 
the next slice since R1 has been automatically updated to 
its proper horizontal position for setting the next bit. The bit 
offset in R1 need only be incremented by + 1 or — 1 to point 
to the exact position of the next bit to be set. 
Diagonal lines are drawn when the value contained in R3 is 
an increment of the bit map's warp. 



example: 

(Assume that the bit map is a 904 x 904 pixel grid.) 



before 



Register Contents 
after 



RO = 1000 RO = 1000 

R1 = 235 R1 = 235 + (150*905) = 135,985 

R2 = 150 R2 = 

R3 = +905 R3 = +905 

This example draws a diagonal line with positive slope start- 
ing at bit position 3 in byte 1029. Notice that the new value 
of R1 = 135,985 is exactly 150 pixels offset from the value 
of R1 in the vertical line drawn in the previous example. 
Adding + 1 to the warp in this example caused the bit posi- 
tion to move not only in the positive vertical direction, but 
also in the positive horizontal direction, forming a diagonal 
line. 

3.2 Implementation of DRAW LINE and SLICE on the 

NS32CG16 

Both a C version of the DRAW LINE algorithm and an 

NS32CG1 6 assembly version are given in the appendix. The 
C program was implemented on SYS32/20 which uses the 
NS32032 processor. An emulation package developed by 
the Electronic Imaging Group at National was used to emu- 
late the SBITS and SBITPS instructions in G, and also the 
MOVMP instruction used for lines longer than 200 pixels. 
The emulation routines, which cover all NS32CG16 instruc- 
tions not available on other Series 32000 processors, are 
available as both G functions and Series 32000 assembly 
subroutines. 

The DRAW LINE program was first written in G using the 

emulation functions. Once this version was tested and func- 
tional, it was translated into Series 32000 code and further 
optimized for speed. The assembly version uses the Series 
32000 assembly subroutines which emulate the SBITS and 
SBITPS instructions. NS32GG16 executable code was de- 
veloped by replacing the emulation subroutine calls with the 
actual NS32GG16 instruction. The functional and optimized 
code was finally executed on the NS32GG16 processor with 
the aid of the DBG1 6 debugger for downloading the code to 
an NS32CG16 evaluation board. Timing for lines of various 
slopes is given in the Timing Appendix. 
Most of the optimization efforts are concentrated in the 
main loop of the SLIGE algorithm. Since the use of SBITS or 
SBITPS for the run length depends on the slope of the line, 
the code is unrolled for the different octants. This minimizes 
branching within the main loop, and cuts down on overall 

execution time. Also, the DRAW LINE takes advantage of 

the NS32CG16's ability to draw fast horizontal, vertical and 
diagonal lines by separating these lines out from the actual 
Bresenham SLICE algorithm. Therefore, time is not wasted 
for trivial lines on executing the initialization sections and 
main loop sections of the SLICE algorithm. 
Branching within the initialization section is also minimized 
by unrolling the code for each octant. Recall from section 
2.2 that in order to extend the algorithm over all octants, the 
abscissa and ordinate displacements must be normalized to 
the first octant and the run length directions must be modi- 
fied to preserve the slope of the line. Partitioning the pro- 
gram into "octant" modules makes the initialization for each 



octant less cluttered with compare-and-branches. Table I 
shows that each octant has a unique value for DELA and 
DELB (the normalized abscissa and ordinate displace- 
ments). Note that at the beginning of the programs, DELX or 
Xf — Xg is checked for sign, and if negative, the absolute 
value function is performed and the starting and ending 
points are exchanged. This is done because each octant 
module of the SLIGE algorithm only cares about the sign of 
DELY with respect to coordinate (Xs,ys). DELX is only impor- 
tant when initializing DELA or DELB, and in this case, only 
the absolute value is needed. 

4.0 SYSTEM SET-UP 

NS32CG16 Evaluation Board: 

— NS32CG16 with a 30 MHz Clock 

—256KB Static RAM Memory (No Wait States) 

— 2 Serial ports 

— M0NCG16 Monitor 

Host System: 

— SYS32/20 running Unix System V 

— DBG16 Debugger 

Software for Benchmarking: 

— START.C Starts timer and calls DRIVER. 

— DRIVER.C Feeds vectors to DRAW_LINE. 

— DRAW LINE.S Line drawing routine which includes 

SLIGE. 
— BIGSET.S Uses MOVMPi to set longer lines. 

Galled by DRAW_LINE if length > 

200. 

4.1 Timing 

Timing Assumptions: 

1 . No wait states are used in the memory. 

2. No screen refresh is performed. 

3. The overhead referred to as the "driver" overhead is the 
time it takes to create the endpoints for each vector. This 
is application dependent, and is not included in the 
Vector/Sec and Pixel/Sec times. 

4. The overhead referred to as the "line drawing" overhead 
is the time it takes to set up the registers for the actual 
line drawing routine. This overhead comes from the 
DRAW LINE program only and is included in all times. 

5. Raw data given in the Timing Appendix for the SBITS, 
SBITPS and MOVMP is the peak performance for these 
instructions. These times do not include line drawing 
overhead or driver overhead. 

The timing for this line-drawing application was done so as 
to give meaningful results for a real graphics application and 
to allow the reader to calculate additional times if desired. 
The routines are not optimized for any particular application. 
All line drawing overhead, such as set-up and branching, is 
included in the given times for Timing Table A, B and G. The 
23 fis driver overhead of the calling routines is not included 
in the given times for vectors per second and pixels per 
second. Calculation of these values was done by subtract- 
ing the 23 fis out of the average time per vector so that the 
given times are only for the processing of the vectors. They 
do not include the overhead of DRIVER.C and START.C 
(refer to these programs in the appendix). 

In addition, the DRAW LINE algorithm is timed for several 

test vectors at various strategic points in the code so that 



the reader may verify set-up times or calculate other rele- 
vant times. The program DRAW LINE.S in the appendix 

contains marl<ers (e.g., T1 , T2 . . . ) for each point at which a 
particular time was taken. The program was run using a 
driver program (DRIVER. C in the appendix) which consists 
of several loops which pass test vectors to the 

DRAW LINE routine. A "return" instruction was placed at 

the time marker so that the execution time was only mea- 
sured up to that marker. These times are given in the Timing 
Appendix Table E and include total execution time up to 
each of the markers. 

A millisecond interrupt timer on the NS32CG16 evaluation 
board was used to time the execution. For each execution, 
the DRIVER program executed its inner loop over 100 
times, and sometimes over 1000 times, so that an accurate 
reading was obtained from the millisecond timer. The final 
times were divided by this loop count to obtain a "bench- 
mark" time. This benchmark time was divided by the total 
number of lines drawn to obtain an average time per vector. 
The overhead of START.C and DRIVER. C in calling the 

DRAW LINE.S routine was not counted in the average 

time per vector or the average time per pixel calculation. 
Table E of the Timing Appendix gives the timing for each of 
the markers and the conditions under which these times 
were taken. 



5.0 CONCLUSION 

The timing for the DRAW LINE algorithm is a good indica- 
tion of the performance of the NS32CG16 in a real applica- 
tion, something which the datasheet specifications can't al- 
ways show. The timing clearly shows that the NS32CG16 is 
well-suited for line-drawing applications. Using the SBITS, 
SBITPS and the MOVMPi instructions, fast line-drawing is 
achieved for lines of all slopes and lengths. The NS32CG16 
is an ideal processor for taking advantage of the much fast- 
er SLICE algorithm. 

The SLICE algorithm, which calculates run lengths of line 
segments to form a complete rasterized line, is much faster 
than its Bresenham predecessor which calculates the line 
pixel by pixel. The SLICE algorithm always executes the 
main loop at least twice as fast as the original Bresenham 
algorithm, which executes its main loop exactly 
max(|delx|,|dely|| times for each line. 
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Series 32000 Programmer's Reference l\/lanual, 1 988. 



Bresenham's SLICE Algorithm: 

1. INITIALIZE PARAMETERS, MAKE NECESSARY ROTATIONS 

2. OUTPUT INITIAL RUN LENGTH (Hq) IN PROPER OCTANT DIRECTION 

MOVE DIAGONALLY IN APPROPRIATE DIRECTION TO START OF NEXT RUN LENGTH 

3. OUTPUT INTERMEDIATE RUN LENGTHS 

COUNT=COUNT-1 

IF COUNT < GOTO 4. 

IFTESTVAR < H = 0-1 AND TESTVAR = TESTVAR + 2*R 

ELSE H = Q ANDTESTVAR = TESTVAR + 2*R-2*DELB 

OUTPUT RUN LENGTH OF LENGTH H IN PROPER DIRECTION 

MOVE DIAGONALLY IN PROPER DIRECTION 

GOTO 3. 

4. OUTPUT FINAL RUN LENGTH OF LENGTH Hp 

5. END 



INITIALIZED PARAMETERS 

DELA = MAXIMUM OF (|DELX|,|DELY|) 

DELB = MINIMUM OF (|DELA|,DELA-MINIMUM||DELX|,|DELY|i 

Q = LOWER[DELA/DELB] 

R = DELA-DELB*Q 

M = LOWER [0/2] 

N = R (IF Q EVEN) 

N = R + DELB (IF O ODD) 

Ho = M (IF DELY>0OR NOO) 

Ho= M-1 (IFDELY<0 AND N = 0) 

Hp = M (IF DELY<0 OR NOO) 

Hp = M-1 (IF DELY >0 AND N = 0) 

COUNT = DELB 

TESTVARo = N + 2*R-2'DELB (IF DELY>0) 

TESTVARo = N + 2*R-2'DELB-1 (IF DELY<0) 



Graphics Image (2000 x 2000 Pixels), 300 DPI 
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FIGURE 6. Star-Burst Benclimarl( 

This Star-Burst image was done on a 2I( x 2k pixel bit map. Each line is 

2I< pixels in length and passes through the center of the image, bisecting the square. The lines are 

25 pixel units apart, and are drawn using the DRAW LINE.S routine. There are a total of 160 lines. 

The total time for drawing this Star-Burst is 1.0s on 15 MHz NS32CG16. 



TIMING APPENDIX 

A. PEAK RAW PERFORMANCE AT 15 MHz 

Function Rate* 

Horizontal Line (SBITS) 9 l\/IBits/s 

Horizontal Line (MOVMP) 60 MBits/s 

Vertical Line (SBITPS) 440 kBits/s 

*Raw performance does not include any register set-up, branching or otlner software set-up overhead. 

B. TRIVIAL LINES (Using Ik x Ik Bit Map Grid) 





Pixels/Line 


Vectors/Sec 


Pixels/Sec 


Comments" 


Horizontal: 


1000 

100 

10 


13,361 
24,136 
45,687 


13,361,838 

2,413,593 

456,870 


Uses BIGSET.S with MOVMP. 
Uses SBITS only. 
Uses SBITS only. 


Vertical and 
Diagonal: 


1000 

100 

10 


424 

3,975 

24,491 


424,000 
397,460 
244,910 


Uses SBITPS. 



**Pixels/Sec and Vectors/Sec are measured from start of DRAW LINE.S only. The 23.128 ;as driver overhead was not included in these measurements. 

C. ALL LINES (Using the "Star-Burst" Benchmark and the SLICE Algorithm) 



Pix/Vector 


Vectors/Sec 


Pixels/Sec 


Total Time* 


Comments** 


1000 

100 

10 


318 
2,811 

14,549 


318,165 
281,118 
145,490 


0.8s 
0.019s 
0.001s 


250 Lines in Star-Burst 
50 Lines in Star-Burst 
10 Lines in Star-Burst 


Avg. Set-up Time Per Line (Measured from Start of DRAW_LINE Only): 37 fis 



D. ALL LINES (Using Original BRESENHAM Iterative Method with SBIT and the Star-Burst Benchmark) 



Pix/Vector 


Vectors/Sec 


Pixels/Sec 


Total Time* 


Comments** 


1000 

100 

10 


163 

1,568 

11,547 


162,746 
158,332 
127,021 


1.5s 
0.033s 
0.001s 


250 Lines in Star-Burst 
50 Lines in Star-Burst 
10 Lines in Star-Burst 


Avg. Set-up Time Per Line (Measured for Line Drawing Routine Only): 30 p.s 



The Bresenham program used for the above table can be found in the Series 32000® Graphics Application Note 5. 

*Total time is measured from start of execution to finish. It includes all line drawing pre-processing, set-up and branching, and it includes all driver overhead of 
DRIVER.C and START.C. This time is a good indication of the pages per minute for the complete Star-Burst benchmark. Vectors/Sec and Pixels/Sec are 

measured from start of DRAW LINE.S only. The 23.712 jis overhead was not included in these measurements. 

**Star-Burst benchmark draws an equal number of lines in each octant. DRIVER.C creates vectors that form the Star-Burst image, passing these vectors to 
DRAW LINE.S as they are created. The bit map image can then be downloaded to a printer for a hard copy, as in Figure 6. 



TIMING APPENDIX TABLE E 


Measurement 
Point 


Measured 
Time/Vector* 


Test Vector 
Used 


Octant of Test Vector 

(Refer to Figure 5) 
And Length of Vector 


Comments 


T1 


23.128 (j,s 


Any Non-Calculated 


Any Octant, Any Length 


Overhead of entry into DRAW_LINE when 
not calculating endpoints of line. Application 
dependent. 


23.712 


STAR-BURST 


All Octants, 1000 Pixels 


Overhead of entry into DRAW_LINE when 
calculating the STAR-BURST vectors. 
Application dependent. 


T2 


40.056 


(0,0,0,999) 


Vertical, 1000 Pixels/Vector 


Average overhead per vertical line 

to start of line draw instruction (SBITPS). 


T3 


41.780 


(0,999,0,0) 


Vertical, 1000 Pixels/Vector 


Average overhead per vertical line with 
negative slope to start of line draw instruction. 


T4 


40.884 


(0,0,999,0) 


Horizontal, 1000 Pix/Vect 


Average overhead per horizontal line to start 
of line draw instruction. (SBITS and BIGSET). 


43.912 


(999,0,0,0) 


Same 


Same as above with negative delta x value. 


T5 


44.532 


(0,0,999,999) 


Diagonal, 1000 Pix/Vect 


Average overhead per diagonal line to start 
of line draw instruction (SBITPS). 


T6 


45.356 


(0,999,999,0) 


Same 


Same as above for diagonal line with 
negative delta x value. 


T7 


71.164 


(0,0,999,10) 


Octant 1 1000 Pix/Vect 


Average overhead per line to first run length 
slice of the SLICE algorithm for octant 1 . 


T8 


87.476 
75.572 
75.568 


(0,0,999,10) 

(0,0,99,10) 

(0,0,9,2) 


Octant 1 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel 
line through first run length of the SLICE 
algorithm. Dependent on the vector length. 


T9 


100.348^3 

88.444 

88.436 


(0,0,999,10) 

(0,0,99,10) 

(0,0,9,2) 


Octant 1 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1 000, 1 00 and 1 pixel 
line to start of main loop of SLICE algorithm. 
Dependent on the vector length. 


T10 


71.856 


(0,0,9,8) 


Octant2 10 Pix/Vect 


Average overhead per line to first run length. 
Not dependent on vector length. 


T11 


79.632 
80.040 
84.180 


(0,0,999,800) 

(0,0,99,80) 

(0,0,9,8) 


Octant 2 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel line 
through first run length of the SLICE algorithm. 
Dependent on the vector length. 


T12 


89.060 
89.476 
105.376 


(0,0,999,800) 

(0,0,99,80) 

(0,0,9,8) 


Octant 2 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1 000, 1 00 and 1 pixel 
line to start of main loop of SLICE algorithm. 
Dependent on the vector length. 


T13 


73.024 


(500,0,700,999) 


Octant 3 1000 Pix/Vect 


Average overhead per line to first run length. Not 
dependent on the vector length. 


T14 


80.736 
80.872 
80.116 


(500,0,700,999) 

(50,0,70,99) 

(5,0,7,9) 


Octant 3 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel line 
through first run length of the SLICE algorithm. 
Dependent on the vector length. 


T15 


89.888 
90.020 
89.268 


(500,0,700,999) 

(50,0,70,99) 

(5,0,7,9) 


Octant 3 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel line 
to start of main loop of SLICE algorithm. 
Dependent on the vector length. 


T16 


73.712 


(10,0,990,999) 


Octant 4 1000 Pix/Vect 


Average overhead per line to first run length. Not 
dependent on the vector length. 


T17 


137.532 

81.148 

78.256 


(10,0,999,999) 

(10,0,90,99) 

(2,0,8,9) 


Octant 4 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel line 
through first run length of the SLICE algorithm. 
Dependent on the vector length. 


T18 


147.236 

90.856 

87.956 


(10,0,999,999) 

(10,0,90,99) 

(2,0,8,9) 


Octant 4 1000 Pix/Vect 

100 Pix/Vect 

10 Pix/Vect 


Average overhead per 1000, 100 and 10 pixel line 
to start of main loop of SLICE algorithm. 
Dependent on the vector length. 


"Each time was measured from start of benchmark execution to the Tx marlcer in the DRAW LiNE.S program. Thus, the overhead of the caiiing routine to the 

DRAW LINE routine isT1 =23.712 fts for the STAR-BURST benchmarlc. Aii programs used for timing are inciuded in the Appendix. Aii times given above are for a 

1k X lie bit map. 



/* This program draws a line in a defined bit map using Bresenham's */ 

/* SLICE algorithm. */ 

i includG<stdio . h> 

Idefine xbytes 25p 

idefine warp 2j3p^ 

fdefine maxy 1999 

unsigned char bit_map[xbYtes*maxy] ; 

extern unsigned char sbitstab[]; 

draw_l ine (xs , ys , xt , yt) 

int xs,ys,xt,yt; 

{ 

int bit, 1, j ,delx,dely,dela,delb, 
hf ,h,hp, testvar,q,r,m, 
n , count , xinc , y inc ; 



delx=xt-xs; 

dely=yt-ys; 

if (xt-xs<^) { 
xs=xt; 
ys-yt; 

delx=abs(delx) ; 
dely= -dely; 
} 

bit=xs-t-ys*warp; 
if (delx==0) { 

if (dely>-)3) { 

sbitps(bit_map, bit, dely, warp) ; 
return ; 
} 
elsef 

sbitps(bit_map,bit,abs(dely) ,-warp) ; 
return ; 
) 
} 
if (dely==P){ 

sbits(bit_map,bit,delx,sbitstab) ; 
return ; 
) 

if (abs(delx)==abs(aely)) { 
if (delx*dely>=0) { 

sbitps(bit_map,bit,abs(dely) ,warp+l) ; 
return ; 
} 
else { 

sbitps(bit_map,bit,delx,-warp+l) ; 
return ; 
} 
} 
if (abs [delx}>abs(dely) ) { 

if (abs(dely)<(delx-abs(dely) ) ) 

dela=delx; 

delb=abs(dely) ; 

xinc=l; 

if (dely>=p) 

yinc^warp ; 
else 

yinc= -warp; 

q=dela/delb; 
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) 

else) 



r=dela-delb*q; 

m=q/2 ; 

if (q-2*(q/2)--p) 

n=r; 
else 

n=r+delb ; 

if ((dely>=P)| I (ni-p)) 

else 

h^=m-l ; 

if ((dely<P) I I (n!=^)) 
hf=m; 

else 

hf=m-l; 

count^delb; 

if (dely>=p) 

testvar=n+2 *r-2 *dGlb ; 
else 

testvar=n+2 *r-2 *delb- 1 ; 
sbits (bit _raap,bit,hp+l , sbitstab) ; 
bit=bit+h0+yinc+xinc ; 

for ( i=count-l ; i>p ; i — ) { 
if (testvar<p) { 

h=q-l; 

testvar+=2*r; 
) 
else { 

h=q; 

testvar+=2*r-2*dGlb; 

} 

sbits (bit_iiiap, bit, h+1, sbitstab) ; 

bit=bit+h+Yinc+xinc ; 

) 

sbits (bit_map, bit, hf, sbitstab) ; 

return; 

dela=abs(delx) ; 

delb=dela-abs(dely) ; 

xinc=l; 

if (dely>=;3) 

yinc=warp; 
else 

yinc= -warp; 
q=dela/delb; 
r=dela-delb*q; 
m=q/2; 
if {q-2*(q/2)==jS) 

n=r; 
else 

n=r+delb ; 
if (CdelY>=^) I I (n!=^)) 

hp=m; 
else 

hp^m-l; 

if ((dely<^) I I (n:=^)) 

hf=m; 
else 

hf =ni-l ; 
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count=delb; 

if Cdely>=^) 

testvar=n+2*r-2*delb; 
else 

testvar=n+2*r-2*delb-l; 
^bitps(bit_Inap,bit,h^^-l,yinc+l) ; 
bit=bit+h^+h^*yinc+l; 

f or ( i=count-l ; i>^ ; i — ) { 
if (testvar<0) { 
h=q-l; 

testvar-l-=2*r; 
} 
else { 

h-q; 

testvar+=2*r-2*delb; 
) 

sbitp£(bit_map,bit,h+l,yinc+l) ; 
bit=bit+h-t-yinc*h+l ; 
) 

sbitps(bit_map,bit,hf+l,yinc+l) ; 
return; 



if (abs(delx)<Cabs(dely)-absCdelx) ) ) ( 
dela=abs(dely) ; 
delb=abs(delx) ; 
yinc=l; 
If (dely>j3) 

xinc=warp ; 
else 

xinc= -warp; 

q=dela/delb; 

r=dela-delb*q; 

m=q/2; 

if Cq-2*(q/2)==0) 

n=r; 
else 

n=r+delb ; 
if ((dely>=0) I I (n!=0)) 

ti0=m; 
else 

hj3=m-l ; 

if (Cdely<j3) | | (n!=0)) 

hf=m; 
else 

hf-m-l; 
count=delb; 

if (dely>=j3) 

testvar=n+2*r-2*delb; 
else 

testvar=n+2*r-2*delb-l; 
sbitps(bit_map,bit,h0+l,xinc) ; 
bit=bit+yinc+ ( i+h^) *xinc; 
f or ( i=CDunt-l ; x>0 ; i — ) { 

if Ctestvar<0) { 
h=q-l ; 
testvar+=2 *r ; 

) 
else { 
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else{ 



h=qf 

testvar-i-=2*r-2*aelb; 

] 

sbitps (bit_inap,bit,h+l,>finc) ; 

bit=bit+yinc+xinc*(l+h) ; 
} 

sbitps (bit_map, bit, hf+l,xinc) ; 
return; 



dela=abs(dely) ; 
delb=dGla-abs (delx) f 
yinc=l; 
If (dGly>p) 

xinc=warp; 
else 

xinc= -warp; 

q=dela/delb; 

r=dela-delb*q; 

in=q/2; 

if (q-2*(q/2)==j3) 

n=r ; 
else 

n=r+delb; 
if ((dely>=p) I I (nl=p)) 

h^=in ; 
else 

h^=m-l; 

if ((dely<P) I I {n!=^)) 

r.f=r[i; 
else 

r.f=ra-l ; 
count=delb; 

if (dely>=p) 

testvar=n+2*r-2*delb; 
else 

testvar=n+2*r-2*delb-l; 
sbitps (bit_raap, bit, hj3+L,xinc+l) ; 
bit=bit+h^+(l+hp)*xinc; 
forCi=count-l; i>^;i — ) { 

If (testvar<^] { 

h-q-l; 

testvar+=2*r r 
) 
else { 

h=q; 

testvar+=2*r-2*delb; 
} 

sbitps (bit_inap,bit,h+l,xinc-i-L) ; 
tlt=bit+h+xinc*{l+h) ; 

sbitps(bit_niap,bit,hf ,xinc+l) ; 
return ; 
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# 


National semiconductor corporation. 




§ 


CTP version 2.4 — draw line.s — Tue Nov 17 13:2p:24 1987 




# 
# 

.file 


compilation options: -0 -S -KC332 -KFj381 -KB4 




"draw_line.s" 






. comm 


bit map,49975p 






.set 


" WARP,2^^^ 






.globl 


draw line 






.globl 


sbitstab 






. align 


4 




_draw_ 


line: 








enter 


[r3,r4,r5,r5,r-7] ,12 




# TL 










movd 


I6(fp),r4 # xf 






movd 


8(fp) ,r5 # xs 






subd 


r5,r4 # delx 






movd 


2^(ep),r6 # yf 






movd 


12(Ep),r7 # ys 






subd 


r7,rG # dely 






cmpqd 


$CP) ,r4 # 0>delx 






ble 


.VERT 






movd 


16(fp),r5 # xf-new xs 






movd 


2^(Ep),r7 # yf=new ys 






absd 


r4,r4 # delx-|delx| 






negd 


r6,r6 # dely=(-dely) 




-VERT 










movd 


r7,ri # ys 






muld 


SWARP,rl # ys*warp 






addd 


r5,rl # bit-ys*WARP+xs 






cmpqd 


$(p) ,ri § delx=p? 






bne 


.HORZ 






cmpqd 


$(p) ,r6 # dely>p? 






bgt 


.VN3G # if no then warp is neg 






addr 


bit map,r^ # set registers for sbitps 






movd 


r6,r2 # r2=dely=length of line 






movd 


$WARP,r3 # r3=warp 




# T2 










sbitps 


# draw line 






exit 


[r3,r4,r5,r5,r7] 
UP) 






ret 






. align 


4 




.VNEG: 










addr 


bi- map,r^ # set reg's for sbitps 






movd 


r6,r2 # r2-C-dely) 






absd 


r2,r2 # r2=dely=length of line 






movd 


$(-WARP),r3 # r3=warp 




# T3 










sbitps 


# draw line 






exit 


[r3,r4,r5,r5,r7] 






ret 






.align 


4 




. HORZ : 










cmpqd 


$(P) ,r6 # aely=p? 






bne 


.DIAG 






addr 


bit niap.r^ # set reg's for sbits 






movd 


r4,r2 # r4=delx=length 






addr 


_sbitstab,r3 # table pointer 




# T4 










sbits 


# try sbits 






bfc 


ok # if not more than 25, skip it 






cmpd 


$2pp,r2 






bit 


bigsl 






addr 


25, r2 






.align 


4 




alpl: 


sbits 








addd 


r2,rl 
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subd 


r2,r4 








cmpd 


r2,r4 








bit 


alpl 








.align 


4 








movd 


r4,r2 








sblts 










exit 


[r3,r4,r5,r6,r7] 








ret 


5(?) 






bigsl: 


bsr 


bigset 






ok: 


exit 


rr3,r4,r5,r6,r7] 








ret 








.align 


4 






. DIAG : 












absd 


r6 , r5 # 


r5=Jdely| 
Idely |=delx? 






cmpd 


r5 , r4 # 






bne 


. SLOPELTl 








cmpqd 


$(?),r6 # 


delY>p? 






bgt 


.DNEG 








addr 


bit inap,rp # 


set reg's for sbitps 






movd 


r4 , r2 # 


r2=delx=length 






movd 


$WARP + l,r3 t 


r3=warp+l for diag 




» T5 












sbitps 


t 


draw line 






exit 


[r3,r4,r5,r6,r7] 








ret 


im 








. align 


4 






.DNEG: 












addr 


bit map,rp # 


set reg's for sbitps 






movd 


r4,r2 # 


r2=delx=lenght 






movd 


$-WARP + l,r3 # 


r3=warp-l for neg slope 




f T6 












sbitps 


* 


draw line 






exit 


[r3,r4,r5,r6,r7] 








ret 


$(*» 








. align 


4 






.SLOPELTl: 


# 


slope less than 1 






cmpd 


r5,r4 # 


Idely |>delx? 






bgt 


.SLOPEGTl 








movd 


r4 , r2 1 


r2=delx 






subd 


r5,r2 t 


delx-|daly| 






cmpd 


r5 , r2 ( 


[dely| >delx-| dely| ? 






bgt 


. 0CTflNT2 1 


If no, start octanti else octant2 






cmpqd 


$(^),r6 t 


dely>.0? 






bgt 


. NEGWARP 








addr 


WARP,-4(fp) I 


pos slope then warp=positive 






br 


.INITl 








.align 


4 






. NEGWARP : 










addr 


-WARP,-4(fp) # 


warp=negative for neg slope 




■INITl 




# 


calculate parameters 






movd 


r4,r3 # 


delx=dela [dely|-delb 






quow 


r5,r3 # 


dela/delb^q 






movd 


r3,r^ # 


calc m 






ashd 


5-1, r^ # 


m=q/2 






movd 


r3,r2 # 


calc r 






mulw 


r5,r2 # 


delb*q 






subd 


r2,r4 # 


r=dela-delb*q 






movd 


r4,r2 # 


set r2 = r 






tbitb 


$^,r3 # 


is r3 odd? 






bfc 


.INIT2 # 


yes, n = r 






addd 


r5,r2 # 


n-r+delb 






.align 


4 






.INIT2 












movd 


r2,r7 # 


pop n 






movd 


r3,tos 1 


push q on stack 






movd 


r^,r2 1 


r2=m=h^ 
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movd 


r^,-8(fp) 


# 


mem=m=hpartb 




cmpqd 


S(?),r7 


# 


n=JI? 




bne 


.INIT3 








cmpqd 


$(?),r6 


# 


dely>^? 




bit 


.INIT4 








addqd 


S-l,r2 


# 


h^=m-l 




br 


.INIT3 








.INIT4: 










subd 


$l,-8(fp) 


# 


hpartb=m-l 




.INIT3: 










addqd 


$l,r2 


# 


takes care of dashes 




addr 


_bit_map,rp 


# 


set reg's for sbits 




addr 


sbitstab,r3 


# 


h^-r2 blt-rl 




1 T7 










sblts 










bfc 


.2 DONE 


# 


set bits if less than 25 




cmpd 


$2W,r2 








bit 


BIGSETl 








movd 


r5, tos 








movd 


r2,r5 








movd 


S25,r2 








.2DRAW25: 










subd 


r2,r5 








sbits 










addd 


r2,rl 








cmpd 


r2,r5 








bit 


.2DR&W25 








movd 


r5,r2 








movd 


tos, r5 








sbits 










br 


.2 DONE 








BIGSETl: 










bsr 


bigsGt 








.2 DONE: 










* T8 










addd 


r2,rl 




blt=bit+h;a+l 




addd 


-4(fp),rl 




bit=bit+h;3+l+warp 




addd 


r4,r4 




2*r 




movd 


r5,r3 




save delb 




addd 


r5,r5 




delb«2 




addd 


r4,r7 




n=n+2*r 




subd 


r5,r7 




testvar=n+2*r+delb*2 




cmpqd 


?(>>), r6 




dely>^ 




bit 


.INIT5 








addqd 


S-l,r7 




testvar-1 




.INIT5: 










movd 


tos,r2 




r2=q=h=run length 
smoothes out line 




addqd 


Sl,r2 






movd 


r3,tos 




push delb^count 




addr 


_sbitstab,r3 




set reg's for sbits 




addr 


bit map,rp 








movd 


-4(£p) ,r6 




warp 




addqd 


S-l,tos 




count=count-l 




cmpqd 


S*l,0(sp) 




count-^? 




bge 


. LASTRUN 








.KSINLOOP: 






Bresenham slice algorithm 




t T9 










cmpqd 


S(»),r7 




testvar>0? 




ble 


. CASE2 








addqd 


$-l,r2 




h=q-l 




addd 


r4,r7 




teEtvar=testvar+2*r 




sblts 










bfo 


.3DRAWLAET 


J 


set bits if less than 25 




cmpd 


$2W,r2 








bit 


BIGEET3 








movd 


r2,tos 
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movd 


r5 , tos 


movd 


r2,r5 


movd 


$25, r2 


.3DRAW25: 




subd 


r2,r5 


sblts 




addd 


r2,rl 


cmpd 


r2,r5 


bit 


.3DRAW25 


movd 


r5,r2 


sbits 




addd 


r2,rl 


movd 


tos,r5 


movd 


tos,r2 


far 


.3 DONE 


BIGSET3 : 




bsr 


bigset 


.3DRAWLAST: 




addd 


r2,rl 


.3 done: 




addd 


r6,rl 


addd 


$l,r2 


addqd 


$(-1) ,tos 


cmpqd 


$m ,?(sp) 


bit 


. MAIKLOOP 


.align 


4 


. LASTRUN : 




cmpqd 


SO) ,tos 


movd 


-8(fp),r2 


sbits 




bfo 


.4 DONE 


cmpd 


S2W,r2 


bit 


BIGSET4 


movd 


r2,tDS 


movd 


r5, toB 


movd 


r2,r5 


movd 


$25, r2 


.4DRRW25: 




subd 


r2,r5 


sbits 




addd 


r2,rl 


cmpd 


r2,r5 


bit 


.4DRAW25 


movd 


r5,r2 


sbits 




addd 


r2,rl 


movd 


tos,r5 


movd 


tos,r2 


br 


. 4 DONE 


BIGSET4 : 




bsr 


bigset 


. 4D0NE : 




exit 


[r3,r4,r5,r6 
$(^) 


ret 


.align 


4 


. CASE2 : 




addd 


r4,r7 


subd 


r5,r7 


sbits 




bfc 


. 5DRAWLAST 


cmpd 


$2W,r2 


bit 


EIGSET5 


movd 


r2,tos 


movd 


r5,tos 


movd 


r2,r5 


movd 


$25, r2 



# update bit 



bit=bit+warp+h+l 
exit h 

count=count-l 
count=^? 



# pop stack 

# hpartb=last run length 

# set bits if less than 25 



# testvar=testvar+2 *r 

# testvar=testvar+2*r-2*dGlb 



# SET BITS IF LESS THAN 25 
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.5DRAW25: 










subd 


r2,r5 








sbits 










addd 


r2,rl 








cmpd 


r2,r5 








bit 


.5DHAW25 








movd 


r5,r2 








sbits 










addd 


r2,rl 








movd 


tos,r5 








movd 


tos,r2 








br 


.5 DONE 








BIGSET5: 










bsr 


bigset 








.5DRAWLAST: 










addd 


r2,rl 


# 


update bit 




.5 DONE: 










addd 


r6,rl 


# 


bit=bit+warp+h4-l 




addqd 


$(-l),tos 


# 


update count 




cmpqd 


$(/3),j3(sp) 


# 


count=p? 




bit 


.MAIKLOOP 








cmpqd 


S(^),tos 


# 


pop stack 




movd 


-8(fp),r2 


# 


hpartb=last run length 




sbits 










bfc 


. 6D0NE 


# 


set bits ir less than 25 




bsr 


bigset 








.6 DONE: 










exit 


[r3,r4,r5,r6 


rV] 






ret 


$(^) 








. align 


4 








. 0CTANT2 : 




f 


draw line in octant 2 




cmpqd 


$(P),r6 


# 


dely>p? 




bgt 


. 2NEGWARP 








addr 


WABP,-4(£p) 


# 


pos slope then warp^positive 




br 


.2INIT1 








.2NEGWARP: 










addr 


-WAKP,-4(£p) 


# 


warp=negative for neg slope 




.2INIT1: 




# 


calculate parameters 




movd 


r4,r3 


t 


dela=delx 




movd 


r2,r5 


« 


delb=delx-|dely| 




quQW 


r5,r3 


# 


dela/delb=q 




movd 


r3,rp 


« 


calc m 




ashd 


$-l,r;! 


f 


m=q/2 




movd 


r3,r2 


# 


calc r 




mulw 


r5,r2 


# 


delb*q 




subd 


r2,r4 


# 


r=dela-delb*q 




movd 


r4,r2 


# 


push r on stack 




tbitb 


$^,r3 








bfc 


.2IN1T2 


# 


then n=r 




addd 


r5,r2 


# 


n=r+delb 




.align 


4 








.2IKIT2; 










movd 


r2,r7 


# 


pop n 




movd 


r3 , tos 


# 


push q on stack 




movd 


r^,r2 


# 


r2=in=h^ 




addqd 


$l,r2 


# 


set one extra bit for smoothness 




movd 


rjl,-8(fp) 


# 


mem=in=hpartb 




cmpqd 


$(^),r7 


# 


n-^? 




bne 


.2INIT3 








cmpqd 


$(^),r6 


# 


dely>^? 




bit 


.2INIT4 








subd 


$l,r2 


# 


hja=m-l 




br 


.2INIT3 








.2INIT4: 










subd 


$1,-8 (fp) 


* 


hpartb=m-l 




.2INIT3: 
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addr 


bit map,r^ 


* 


set reg's for sbits 




movd 


-4(fp),r3 


f 


warp-r3 hp=r2 bit=rl 




addqd 


$l,r3 


* 


octant 2 needs diag runs 




* TljS 










sbltps 




# 


draw first run length 




» Til 










addqd 


$l,rl 


# 


update bit in x direction 




subd 


r3,rl 


# 


sbitps adds extra warp 




addd 


r4,r4 


* 


2*r 




movd 


tos , r2 


« 


q=h=next run length 




addqd 


$l,r2 


t 


set extra bit for smoothness 




movd 


r5 , tos 


I 


push delb=count 




addd 


r5,r5 


f 


delb*2 




addd 


r4,r7 


t 


n=n+2*r 




subd 


r5,r7 


# 


testvar-n+2 *r+delb*2 




cmpqd 


$(*)),r6 


# 


dely>^ 




bit 


.2INIT5 








subd 


$l,r7 


# 


testvar-1 




.2INIT5: 










subd 


$l,tos 


« 


count=count-l 




cmpqd 


??,?(sp) 


# 


count^p? 




bge 


.2LASTRUN 








.2MAINL00P: 




« 


Bresenham slice alqorithm 




# T12 










cmpqd 


$(«>), r? 


# 


testvar>^? 




ble 


. 2CASE2 








subd 


$l,r2 


» 


h=q-l 




addd 


r4,r7 


# 


testvar=testvar+2*r 




movd 


r2,tos 


# 


preserve h 




sbitps 




# 


draw diag line of length h 




movd 


tos,r2 


# 


renew h 




addqd 


$l,rl 


# 


update bit in x direction 




subd 


r3,rl 


# 


sbitps adds one warp extra 




addd 


Sl,r2 


# 


exit h to q 




subd 


SI, tos 


# 


count=count-l 




cmpqd 


S?,?(sp) 


# 


count^^? 




bit 


. 2MAINL00P 








.align 


4 








.2LASTRUN: 










cmpqd 


$(P),tos 


# 


pop stack 




movd 


-8(fp),r2 


# 


hpartb=last run length 




sbitps 




# 


all other reg's set up 




exit 


[r3,r4,r5,r6 


r7] 






ret 


$(*» 








-align 


i 








.2CASE2: 










addd 


r4,r7 


# 


testvar=testvar+2 *r 




subd 


r5,r7 


# 


testvar=testvar+2*r-2*delb 




movd 


r2,tos 


# 


preserve h 




sbitps 




# 


draw line of length h=q 




movd 


tos,r2 


# 


renew h 




addqd 


$l,rl 


# 


update bit in x direction 




subd 


r3,rl 


# 


sbitps adds one warp extra 




subd 


S(l) ,tos 


# 


update count 




cmpqd 


S^.;l(sp) 


# 


count=^? 




bit 


. 2MAINU)0P 








cmpqd 


S(^) ,tos 


# 


pop stack 




movd 


-8tfp),r2 


# 


hpartb=last run length 




sbitps 




# 


all other reg ' s set up 




exit 


rr3,r4,r5,r6 


r7] 






ret 








.align 4 








.SLOPEGTl: 




# 


coordinates are rotated for these lines 




movd 


r5,r2 


# 


r2=Jdely| 
1 deiy 1 -delx 
delx> 1 dely | -delx? 




subd 


r4,r2 


# 




cmpd 


r4,r2 


# 
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bgt 


.20CTANT2 


# 


if no, start octantl else octant2 




cmpqd 


S(P),r6 


# 


delY>^? 




bgt 


. 3NEGWAEP 








addr 


WAEP,-4(fp) 


# 


pos slope then warp=positive 




br 


.3INIT1 








. 3NEGWARP : 










addr 


-WAEP,-4(fp) 


# 


warp=nGgative for neg slope 




.3INIT1: 




t 


calculate rotated parameters 




movd 


r5,r3 


* 


dela=|dalY| 
delb=delx 




movd 


r4,r5 


# 




movd 


r3,r4 


* 


dela in r4 




quow 


r5,r3 


# 


dela/delb=q 




movd 


r3,r^ 


* 


calc m 




ashd 


S-I,rjl 


« 


m=q/2 




movd 


r3,r2 


« 


calc r 




mulw 


r5,r2 


« 


delb*q 




subd 


r2,r4 


t 


r=dela-delb*q 




movd 


r4,r2 


« 


push r on stack 




tbitb 


S^,r3 








bfc 


.3INIT2 


# 


then n=r 




addd 


r5,r2 


# 


n=r+delb 




.align 


4 








.3INIT2: 










movd 


r2,r7 


# 


pop n 




movd 


r3,tos 


# 


push q on stack 




movd 


rp,r2 


* 


r2-m.h^ 




addqd 


Sl,r2 


# 


set one extra bit for smoothness 




movd 


r^,-S(fp) 


# 


mem=m=hpartb 




cmpqd 


S(?),r7 


# 


n=^? 




bne 


.3INIT3 








cmpqd 


S(^),l-6 


# 


delY>^? 




bit 


.3INIT4 








subd 


$l,r2 


* 


h)J=m-l 




br 


.3INIT3 








.3INIT4; 










subd 


Sl,-8(fp) 


* 


hpartb=m-l 




.3INIT3: 










addr 


bit map.r^ 


« 


set reg's for sbits 




movd 


-4(fp),r3 


« 


warp=r3 hp=r2 bit^rl 




J T13 










sbitps 




* 


draw first run length 




i T14 










addqd 


$l,rl 


# 


update bit in x direction 




addd 


r4,r4 


# 


2*r 




movd 


tos,r2 


# 


q=h=next run length 




addqd 


$l,r2 


# 


set extra bit for smoothness 




movd 


r5,tos 


* 


push delb=count 




addd 


r5,r5 


* 


delb*2 




addd 


r4,r7 


# 


n=n+2*r 




subd 


r5,r7 


# 


testvar=n+2*r+delb*2 




cmpqd 


S(^),r6 


t 


dely>p 




bit 


.3INIT5 








subd 


Sl,r7 


f 


testvar-1 




.3INIT5: 










subd 


Sl.tos 


* 


count=count- 1 




cmpqd 


iH.Hisp) 


i 


count-^? 




bge 


.3LASTRUN 








.3MAINL00P: 




# 


Bresenham slice algorithm 




* T15 










cmpqd 


5(?),r7 


# 


testvar>^? 




ble 


.3CASE2 








subd 


$l,r2 


t 


h=q-l 




addd 


r4,r7 


t 


testvar=testvar+2*r 




movd 


r2,tos 


t 


preserve h 




sbitps 




t 


draw vert line of length h 




movd 


tos,r2 


# 


renew h 
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addqd 


$l,rl 


# 


update bit in x direction 




addd 


$l,r2 


1 


exit h to q 




subd 


Sl,tos 


# 


count=count-l 




cmpqd 


S^,?(sp) 


t 


count=^? 




bit 


. 3MAINL00P 








•align 


4 








.3LASTRUN: 










cmpqd 


SC*1) ,tos 


# 


pop stack 




itiovd 


-8(fp),r2 


# 


hpartb=last run length 




sbitps 




# 


all other reg's set up 




exit 


[r3,r4,r5,r6 


r7] 






ret 








.align 


4 








.3CASE2: 










addd 


r4,r7 


# 


testvar=testvar+2 *r 




subd 


r5,r7 


# 


testvar=testvar+2*r-2*delb 




movd 


r2,tos 


# 


preserve h 




sbitps 




# 


draw line of length h=q 




movd 


tos,r2 


♦ 


renew h 




addqd 


$l,rl 


# 


update bit in x direction 




subd 


$(1) ,tos 


# 


update count 




cmpqd 


5^,^(sp) 


# 


count=p? 




bit 


. 3MAINL00P 








cmpqd 


$(^),tos 


# 


pop stack 




movd 


-8(fp),r2 


# 


hpartb=last run length 




sbitps 




# 


all other reg's set up 




exit 


[r3,r4,r5,r6 


r7] 






ret 








.align 


4 








. 20CTANT2 : 




# 


draw line in octant 2 




cmpqd 


H?],r6 


# 


dely>^? 




bgt 


. 4NEGWARP 








addr 


MARP,-4 [fp) 


# 


pos slope then warp=positive 




br 


.4INIT1 








.4NEGHARP: 










addr 


-WARP, -4 (fp) 


* 


warp-negative for neg slope 




.4INIT1: 




# 


calculate parameters 




movd 


r5,r3 


1 


dela=delx 




movd 


r5,r4 


f 


dela into r4 




movd 


r2,r5 


* 


delb-delx-ldelyl 




quow 


r5,r3 


# 


dela/delb=q 




movd 


r3,r0 


1 


calc m 




ashd 


?(-!) ,rp 


* 


m=q/2 




movd 


r3,r2 


# 


calc r 




mulw 


r5,r2 


* 


delb»q 




subd 


r2,r4 


# 


r=dela-delb*q 




movd 


r4,r2 


1 


push r on stack 




tbitb 


??.r3 








bfo 


.4INIT2 


I 


then n=r 




addd 


r5,r2 


# 


n=r+delb 




.align 


4 








.4INIT2: 










movd 


r2,r7 


1 


pop n 




movd 


r3,tos 


# 


push q on stack 




movd 


r?,r2 


1 


r2=m=h^ 




addqd 


$l,r2 


* 


set one extra bit for smoothness 




movd 


r;i,-8(fp) 


# 


inem=in=hpartb 




cmpqd 


$(«»,r7 


1 


n=?? 




bne 


.4IKIT3 








cmpqd 


$(?) ,r6 


# 


dely>JI? 




bit 


.4INIT4 








subd 


$l,r2 


# 


h^=in-l 




br 


.4INIT3 








.4INIT4: 










subd 


51,-8(fp) 


♦ 


hpartb=in-i 




.4INIT3: 
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addr 


bit map , r^ 


# 


set reg's for sbits 




movd 


-4(fp),r3 


# 


warp-r3 hp=r2 bit-rl 




addqd 


$l,r3 


# 


octant 2 needs diag runs 




# T16 










sbitps 




# 


draw first run length 




# T17 










subd 


$l,rl 


# 


update bit 




addd 


r4,r4 


# 


2*r 




movd 


tos,r2 


# 


q=h=next run length 




addqd 


$l,r2 


# 


set extra bit for smoothness 




movd 


r5,tos 


# 


push delb=count 




addd 


r5,r5 


# 


delb*2 




addd 


r4,r7 


# 


n=n+2*r 




subd 


r5,r7 


# 


testvar=n+2*r+delb*2 




cmpqd 


S(^),r6 


# 


dely>0 




bit 


.4INIT5 








subd 


51, r? 


# 


testvar-1 




.4INIT5: 










subd 


$l,tos 


# 


count=count-l 




cmpqd 


S^,^(sp) 


# 


count =^? 




bge 


. 4LASTRUN 








.4MAINL00P: 




# 


Bresenham slice algorithm 




# TIS 










cmpqd 


S»<),I-7 


# 


testvar>p? 




ble 


. 4CASE2 








subd 


Sl,r2 


# 


h=q-l 




addd 


r4,r7 


# 


testvar=testvar+2*r 




movd 


r2 , tos 


# 


preserve h 




sbitps 




# 


draw diag line of length h 




movd 


tos,r2 


# 


renew h 




subd 


$l,rl 


# 


sbitps adds one warp extra 




addd 


$l,r2 


# 


exit h to g 




subd 


SI, tos 


# 


count=count-l 




cmpqd 


5;i,p(sp) 


# 


count-p? 




bit 


.4HAINLO0P 








.align 


4 








.4LASTRUN: 










cmpqd 


$iP),tos 


# 


pop stack 




movd 


-8(fp),r2 


# 


hpartb=last run length 




addqd 


$l,r2 








sbitps 




# 


all other reg's set up 




exit 


[r3,r4,r5,r5 


r7] 






ret 








.align 


4 








.4CASE2: 










addd 


r4,r7 


# 


testvar=tGstvar+2 *r 




subd 


r5,r7 


# 


testvar=testvar+2*r-2*delb 




movd 


r2,tos 


# 


preserve h 




sbitps 




# 


draw line of length h-q 




movd 


tos,r2 


# 


renew h 




subd 


$l,rl 


# 


sbitps adds one warp extra 




subd 


$(l),tOB 


# 


update count 




cmpqd 


S?,«l(sp) 


# 


count=^? 




bit 


.4HAINL00P 








.align 


4 








cmpqd 


$(?) ,tos 


♦ 


pop stack 




movd 


-8(fp),r2 


# 


hpartb=last run length 




addqd 


$l,r2 








sbitps 




# 


all other reg ' s set up 




exit 


rr3,r4,r5,r6 


r7] 






ret 
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.globl 


bigset 


save 


[r^,rl,r2,r3,r4,r5,r6] 


movd 


rl,r4 


ashd 


S-3,r4 


addcJ 


r4,r? 


andd 


57, rl 



BIGSET. S uses MOVMP and the OR instructions to set long horizontal lines 

ff 

bigset: save [r^ , rl, r2 , r3 , r4 , r5, r6] #save registers wc will affect 

#get current bit offset 
idivide by eight to get byte offset 
#add in base, rji is new base pointer 
#niask off msb ' s of bit pointer to 
#get bit = bit offset mod 8 

#Now we have true base address and bit offset within base. Now we will move 
#to double word alignment. This speeds up the MOVMPD for long bit sequences. 

#place mask in r4 

#gQt low two bits of address 

#and get bytes left to alignment 

#rem += 1 (for the byte we are on) 

#rem *= 8 to get bits to alignment 

#subtract current bit offset 

#is this more than number of bits left 

#it is, do it the short way 

#if WG are already double aligned, go 

#do the MOVMPD 

^calculate index into table 

#index = 32 * bit offset 

#index += run length 

#or in required bits 

#clear last two bits, and 

#bump to next double 

#zap sp'd bits off 

#save run length for a minute 

#and save pointer to table 

#rl ^ rl / 32 = number of doubles 

#^et source pattern from table 

# increment is rl 

# yes, use instruction 

#mask off all but last 32 bits 

#insert the last few bits 

^restore saved registers 



#check to see if it is exactly 
#32 bits. If it is, branch. 
^calculate index into table 
#index = 32 * bit offset 
#index += run length 
#or in required bits 
#restore saved registers 



movqd 


3,r4 


andd 


r/i),r4 


xorb 


$3,r4 


addqd 


l,r4 


ashd 


$3,r4 


subd 


rl,r4 


cmpd 


r4,r2 


bge 


shrt 


cmpd 


$32, r4 


beq 


Tnvm 


movd 


rl,r5 


Ishd 


$5,r5 


addd 


r4,r5 


ord 


r3[r5:d],^(rj3) 


bleb 


$3,rjl 


addqd 


4,rp 


subd 


r4,r2 


movd 


r2,r4 


movd 


r3,r5 


ashd 


S-5,r2 


movd 


1^2^(r3),r3 


movqd 


4,rl 


movmpd 
andd 


$Jlxlf ,r4 


ord 


r5[r4:d],p(rp) 


restore 


[r^,rl,r2,r3,r4,r5,r61 


ret 


$^ 


. align 


4 


cmpb 


532, r2 


beq 


shrtl 


movd 


rl,r4 


Ishd 


$5,r4 


addd 


r2,r4 


ord 


r3[r4:d],JI(r^) 


restore 


[ rcl , rl , r2 , r3 , r4 , r5 , r6 1 


ret 


S*> 


.align 


4 


movd 


lj)2^(r3),^(r^) 


restore 


r r^ , rl , r2 , r3 , r4 , r5 , r6 ] 
5^ 


ret 



#copy last entry of table 
#Call 32 bits) and restore 
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0) 

o 



(0 



0) 



/* Program driver. c feeds line vectors to LINE_DRAW.S forming Star-Burst. 



#lnclude <stdio.h> 
#define xbytes 2 5^ 
#definQ maxx 1999 
(/) #define maxy 1999 



Q. 

i0 for (count=l;count<=l^^^;test++) { 



unsigned char bit_map[xbytes*maxy] ; 

mainO 

{ 

int i , count ; 

/* generate Star-Burst image */ 



for (i=P;i<=iiiaxY;i+=25) 

draw_line(p, i,ma 
for {i=^;i<=maxx;i+=25) 

draw line (i, maxy ,itiaxx-i ,^) ; 



Or or I i=p; i-^^maxy; i-^=zoj 
draw_line(p, i,maxx,inaxy-i) 

O 

CM 
CO 
(/) 



(^ /* Start timer and call main procedure of DRIVER. C to draw lines */ 

O start { 
long *timer ^ (long *) px6^p ; 
- k *timGr ~ ^; /* write a zero to timer location */ 

^J mainC^,^); /* Show argc as zero, argv ->^ */ 

^^ return ( *timer) ; /* return, in rp, the current time */ 

CO 

(/) 

Z 
0) 
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LIFE SUPPORT POLICY 

NATIONAL'S PRODUCTS ARE NOT AUTHORIZED FOR USE AS CRITICAL COMPONENTS IN LIFE SUPPORT 
DEVICES OR SYSTEMS WITHOUT THE EXPRESS WRITTEN APPROVAL OF THE PRESIDENT OF NATIONAL 
SEMICONDUCTOR CORPORATION. As used herein: 

1. Life support devices or systems are devices or 2. A critical component Is any component of a life 

systems which, (a) are Intended for surgical Implant support device or system whose failure to perform can 

into the body, or (b) support or sustain life, and whose be reasonably expected to cause the failure of the life 

failure to perform, when properly used In accordance support device or system, or to affect Its safety or 

with Instructions for use provided in the labeling, can effectiveness, 
be reasonably expected to result In a significant Injury 
to the user. 

^\ National Semiconductor National Semiconductor National Semiconductor National Semiconductor 

^U ^^m Corporation Europe Hong Kong Ltd. Japan Ltd. 

f^ I/' 1 1 1 1 West Bardin Road Fax: ( + 49) 0-180-530 85 86 13th Flodr, Straight Block, Tel: 81-043-299-2309 

^jj Arlington, TX 76017 Email: cnjwge@tevm2.nsc.com Ocean Centre, 5 Canton Rd. Fax:81-043-299-2408 

lO Tel: 1 (800) 272-9959 Deutsch Tel: ( + 49) 0-180-530 85 85 Tslmshatsul, Kowloon 

Fax: 1(600) 737-7018 English Tel: (+49) 0-180-532 78 32 Hong Kong 

Frangais Tel: ( + 49) 0-180-532 93 58 Tel: (852) 2737-1600 

Italiano Tel: (+49) 0-180-534 16 80 Fax:(852)2736-9960 

National does not assume any responsibility lor use oi any circuitry descnbed, no circuit patent iicenses are implied and National reserves the right at any time without notice to change said circuitry and speciiications. 



